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We calculate the O(a) improvement coefficient for the axial-vector current using the nonperturbative method 
proposed by the LANL group. Results for the coefficient in the range (3 — 5.93 to 6.2 are presented. We find ca is 
close to the 1-loop tadpole-improved perturbative value. In addition, scaling of the pseudoscalar decay constant 
and renormalised quark mass is improved compared to that obtained using the (larger) ca values obtained by the 
ALPHA collaboration. 



1. Introduction 

Symanzik improvement offers reduced discreti- 
sation errors for masses and matrix elements at 
finite lattice spacing but requires the calculation 
of the coefficients which accompany the improve- 
ment terms. In the light quark regime these co- 
efficients can be calculated nonperturbatively by 
imposing the chiral Ward identities on the lattice. 

ca is the 0(a) improvement coefficient for the 
axial- vector current. 



+ acAdi + 0(a 2 ) 



(1) 



It has been calculated nonperturbatively by the 
ALPHA collaboration @ and LANL group @. 
The results are summarized in figure 0. There 
is a significant discrepancy in the values of ca 
at smaller (3. This might be accounted for by 
the 0(a) ambiguity that exists in ca when it 
is determined nonperturbatively. Nevertheless, 
larger values of ca lead to stronger lattice spac- 
ing dependence of and the quark mass, and 
increase the difficulty of performing a continuum 
extrapolation. In addition, the ALPHA results 
have a very different (3 dependence compared to 
the 1-loop perturbative values, which suggests 
they contain a large nonperturbative contribu- 
tion. Our aim was to investigate how well ca 
is determined using the LANL method, which 
only requires conventional analysis, rather than 
Schrodinger functional techniques of the ALPHA 
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collaboration. 

ca is determined nonperturbatively using the 
PCAC relation. In euclidean space, for zero mo- 
mentum 



< d i A i (x)J >= 2mpcAC < P(x)J > 



(2) 



where mpcAC is the PCAC quark mass, P = 
ip"f5tp and J is any operator with the pseudoscalar 
quantum numbers. Defining rj(t) = <d ^p^')^ > , 
equation @ becomes on the lattice 



rj(t) = m(t) 



2m PC AC + O(a). 



(3) 



The size of the 0(a) terms depends on the num- 
ber of states contributing to rj(t). At early 
times (t ex ), the significant contribution from ex- 
cited states means the O(a) term is larger than 
at later times when the ground state domi- 
nates (t gs ). The improvement of equation [l] re- 
moves the 0(a) terms for all states and hence the 
two masses 



rj(t gs ) + ac A sj(t gs ) 

rj{t e x) + ClCASj{t ex ) 



(tgs) 



2m PC AC + 0{a 2 ). (4) 
2m PC AC + 0'{a 2 ).(5) 



where sj(t) — -^^rp^r > are equal up to higher 
order terms. We can use these equations to solve 
for ca- 

1 r(t ex ) -r{t gs ) _ 

7T~\ 77 \ = CA 6 

a s(t ex ) - s(t gs ) 

The LANL method involves performing a fit to 
rj(t)+acASj(t) — constant ~ 2m over the fitting 
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Figure 1. Summary of ca determinations. LANL- 
2 and LANL-3 refer to ca determined using 2- 
point and 3-point lattice derivatives respectively, 
as described in reference . 

range t ex to t gs and is numerically equivalent to 
equation |[ ca and 2m are extracted from the fit. 
The advantage of performing a fit is that one can 
test the ansatz with the x 2 . 

The 0(a) ambiguity in ca is due to 0(a 2 ) terms 
in the axial- vector current, the pseudoscalar cur- 
rent and the light quark action. We define ca in 
the chiral limit, but at finite quark mass there 
is an additional source of O(a) error if standard, 
symmetric, 0(a 2 ) lattice temporal derivatives are 
used. This error can be removed by using im- 
proved 0(a 4 ) lattice derivatives. While, both 
choices of lattice derivative lead to consistent val- 
ues for ca in the chiral limit, we found it advan- 
tageous to use improved derivatives as this weak- 
ened the dependence of ca on the quark mass 
and made the chiral extrapolation easier. We also 
found a wider window in t ex from which ca can 
be extracted. 

2. Results for ca 

We calculated ca on the UKQCD Collabora- 
tion quenched data set, at j3 = 5.93, 6.0 and 6.2; 
simulation details can be found in |3j. The best 
analysis was possible at = 5.93. 

Table [j] presents the results for ca at = 5.93 
for "fuzzed" source || and LL correlators, where 
improved lattice derivatives have been used. We 
see that there is agreement between the results 
from different smearings and the results are stable 
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k = 0.1327 




FL 


4 


-0.037(2) 


-0.050(3) 


5 


-0.035(4) 


-0.052(8) 


6 


-0.029(9) 


-0.022(19) 


7 


-0.025(15) 




8 


-0.066(264) 




LL 


6 


-0.024(6) 


-0.032(14) 


7 


-0.022(11) 




8 


-0.046(30) 





Table 1 

ca values extracted for = 5.93 at finite quark 
mass and extrapolated to the chiral limit. For all 
results shown Q for the fit is > 0.01. At k c and 
higher t ex the statistical errors dominate and the 
results are not shown. t gs = 14. 

with t ex . However, there is only a small window of 
timeslices for which Q for the fit is reasonable and 
before noise dominates. No significant change in 
ca was found when t gs was varied. 

We take ca = —0.032(14) as our final value for 
this 0. In choosing a value with a conservative er- 
ror we aim to take into account the residual 0(a) 
ambiguity in ca and some of the other system- 
atic errors. In principle all results in the table are 
valid estimates of ca- However, choosing a result 
with a much smaller error would require much 
more care in keeping the improvement conditions 
applied fixed when ca is recalculated at other /3s; 
thus ensuring, for example, any residual 0(a) er- 
ror in ca vanishes when the continuum limit is 
taken. In particular, with the LANL method one 
would need to keep the relative contribution of 
ground state to excited states at t ex (in physical 
units) constant and this requires tuning of the 
smearing function, which was not possible in this 
study. 

Repeating the analysis at = 6.0 we found 
evidence of significant finite volume dependence. 
ca decreased from —0.039(6) at k = 0.13344 on a 
16 3 lattice to -0.003(10) on a 32 3 lattice. The 
16 3 lattice is 16% smaller in extent than that 
used at (3 — 5.93. Due to problems with chiral 
extrapolations, our final results as a function of 
t ex either had too small a statistical error to re- 
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Figure 2. Scaling of the renormalised PS decay 
constant in units of rp. 
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Figure 3. Scaling of the renormalised quark mass 
in units of ro 



fleet the change in volume, or had too large an 
error to be useful. Hence, we disregard the re- 
sults at this (3 in the later analysis. At (3 = 6.2, 
the lattice size is only 6% smaller, and we obtain 
ca — —0.031(5). We assume the latter error is 
sufficient to compensate for the small change in 
volume from (3 = 6.2 to 5.93. 

3. Scaling of f v and m M 

We calculate the renormalised pseudoscalar de- 
cay constant, given by 

f rm = Z A {\ + ab A m q )(fW + ac A f {1) ) (7) 
/(°) =< 0\Ai\P >, /W =< 0\d±P\P > (8) 

using our results for ca and nonperturbative de- 
terminations of Za and the perturbative results 
for bA', f ren is extracted at a reference mass 
(roMps) 2 — 3.0. This is compared in figure || 
with the results using the ALPHA collaboration 
values for ca- We see that scaling is much better 
when our smaller values of ca are used. The si- 
multaneous fit to all 4 points only has a marginal 
X 2 = 2.4. More points and higher orders in a 
would be needed to include the ALPHA ca val- 
ues for f ren . 

The renormalised quark mass, at the same ref- 
erence mass, was calculated using the method of 
reference The renormalisation-group invari- 
ant mass is computed from the PCAC mass us- 
ing the nonperturbatively determined renormal- 
isation factor Zm'- M = Zm^pcac The MS 
mass at 2 GcV is then obtained using 4-loop per- 
turbation theory. Figure || shows that the scaling 



is better for smaller values of ca', the x 2 for the 
combined fit is 3.3. 

4. Conclusions 

We implemented the LANL method for deter- 
mining ca in the range j3 = 5.93 to 6.2. Figure |l| 
shows our results are consistent with those of the 
LANL group. In particular, ca does not change 
rapidly with j3 and is close to the 1-loop tadpole- 
improved value (using ap(l/a)). At {3 = 6.0, 
evidence of a finite volume dependence in ca in- 
dicated care must be taken in keeping the physical 
volume in the simulation fixed when changing (3 
with this method. 

The scaling of both and m MS are improved 
using our, smaller, values of ca compared to the 
values obtained by the ALPHA collaboration. 
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